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Numerical Aperture of 
Single-Mode Photonic Crystal Fibers 

Niels Asger Mortensen, Jacob Riis Folkenberg, Peter M. W. Skovgaard, and Jes Broeng 
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Abstract — We consider the problem of radiation into free 
space from the end-facet of a single-mode photonic crystal 
fiber (PCF). We calculate the numerical aperture NA = sin 6 
from the half-divergence angle 9 ~ tan" 1 (X/ttw) with ttw 2 be- 
ing the effective area of the mode in the PCF. For the fiber 
first presented by Knight et al. we find a numerical aperture 
NA ~ 0.07 which compares to standard fiber technology. We 
also study the effect of different hole sizes and demonstrate 
that the PCF technology provides a large freedom for NA- 
engineering. Comparing to experiments we find good agree- 
ment. 

Keywords — Photonic crystal fiber, numerical aperture, 
Gaussian approximation 



(N 
> 
CO 

o 

O 
CI 

o 

O 



I. Introduction 

PHOTONIC CRYSTAL FIBERS (PCF) constitute a 
completely new class of optical fibers consisting of pure 
silica with air-holes distributed in the cladding. Among 
many remarkable properties |Q PCFs are believed to have 
a potential for high-numerical aperture (NA) applications. 
Here we report a calculation of the NA for the class of PCFs 
first fabricated by Knight et al. [0 , || . For this particular 
fiber we find a numerical aperture up to NA ~ 0.07. We 
also demonstrate how the NA may be controlled by the 
hole size for a given pitch and wavelength. 




Fig. 1 

Coupling of light from end-facet of fiber (z = 0) into free 
space. The insert shows a micrograph of the end-facet of a 

PCF. 



The paper is organized as follows: First we consider the 
problem of radiation into free space from the end-facet of a 
single-mode optical fiber with a mode approximated by a 
Gaussian of width w. Solving the scattering problem at the 
end-facet of the fiber exactly we check the range of validity 
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of the text-book result 8 ~ tan _1 (A/7ru?) (see e.g. Ref. Q) 
where A c s = to 2 is the effective area. For w<A we find 
deviations whereas nice agreement is found for w > A. We 
then turn to the application of the PCF of Knight et al. 
H which belongs to the latter regime with w > A. Finally, 
we compare our calculations to experiments. 

II. Numerical aperture in the Gaussian 

APPROXIMATION 

The numerical aperture NA = sin 9 (see Fig. 0) may 
be defined in various ways, but often one defines it in the 
far-field limit (z — ► oo) from the half-divergence angle 6 V 
between the z-axis and the z^- intensity point r v {z), i.e. 



tan#„ = lim 



r v {z) 



with r v [z) determined from 

\^ > (z,r± = r v ) 



\9 > {z,r J _ = 0)\ 



(1) 



(2) 



For a Gaussian field \t of width w one has the standard 
approximate expression for v = 1/e 2 ~ 13.5% B 



2 A 
tanc' 1 / e 2 ~ - — — . 

kw TTW 



For the v = 5% intensity point 



tan # 5% = 



In 20 



x tan( 



Ve 



(3) 



(4) 



which is often the one used experimentally. Eqs. (|3j,[4j) are 
valid for kw 3> 1, but in order to check the validity in 
the limit with kw of order unity we solve the scattering 
problem at the end-facet of the fiber exactly. In the fiber 
(z < 0) the field is of the form 

#<(r) oc iP(r ± )(e ll3{uj)z +^e-' l ^ z ), z<0 (5) 

where the transverse field is approximated by a Gaussian 

! / ) (r 1 )oce-( r ^" , ' 2 I (6) 

which has an effective area A e g = irw 2 at frequency to. At 
the end-facet of the fiber (z = 0) the field couples to the 
free-space solution 



*>(r) oc J dk ± ^(k ± )( 



z>0 



(7) 
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which is a linear combination of plane waves with lu = ck = 
c(2ir/\) and k = k± + k\\. 

In order to solve the elastic scattering problem, Auj — 
U)(f3)—u)(k) = 0, we apply appropriate boundary conditions 
at the end-facet of the fiber; continuity of VP and d^> /dz. At 
z = we thus get two equations determining the reflection 
amplitude Si and the transmission amplitude S? '. Elimi- 
nating and substituting the resulting ST into Eq. (0) we 
get 



2.5 



*>(r) oc 27rfc 2 



dxx 



2n r 



«cff 



xe -(xkw/2) 2 j a ( xkr ^yV^ k * . (8) 

Here, x = k±/k, Jq is the Bessel function of the first kind 
of order 0, and n e g — (3/k is the effective mode-index. 
Eq. (||) is the exact solution to the scattering problem and 
in contrast to many approximate text-book expressions (see 
e.g. Refs. j|) we have here treated the scattering problem 
correctly including the small, but finite, backscattering in 
the fiber. Thus, we take into account the possible filtering 
in transmitted k± at the fiber end- facet. The solution has 
similarities with the Hankel transform usually employed 
in the far-field inversion integral technique, see e.g. p|. 
Numerically we have found that Eq. (|§|) gives a close-to- 
Gaussian field in the far- field limit. 

In Fig. U we compare the two approximate solutions 
Eqs. (p^) to a numerically exact calculation of tan 9 V from 
Eq. (H). The calculation is performed for the realistic sit- 
uation with n c ff = (3/k = 1.444 corresponding to a silica- 
based fiber. For kw ~ 1 the deviations increase because 
of the small, but finite, backscattering at the end-facet of 
the fiber. For kw somewhat larger than unity a very nice 
agreement is found. A typical all-silica fiber like the Corn- 
ing SMF28 has kw > 10. 

III. Application to photonic crystal fibers 

We consider the class first studied in Ref. || which con- 
sists of pure silica with a cladding with air-holes of diameter 
d arranged in a triangular lattice with pitch A. For a review 
of the operation of this class of PCFs we refer to Ref. || . 

In applying Eq. ^ to PCFs we calculate w from the 



effective area A c s 



■kw given 



by 



[fdr\H(r, 



Jdr\H(r,z)\ 



(9) 



Indeed we find that the corresponding Gaussian of width 
w accounts well for the overall spatial dependence of the 
field. Of course we thereby neglect the satellite spots seen 
in the far-field but because of their low intensity they 
only give a minor contribution to the NA j^] . 

For the field H of the PCF, fully-vectorial eigenmodes 
of Maxwell's equations with periodic boundary conditions 
are computed in a planewave basis || . 

Figure H illustrates the effective mode- index and effective 
area as a function of wavelength for a PCF with d/A = 
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Fig. 2 

Plot of tan as a function of the dimensionless parameter 
kw. The points are the results of a numerical exact 

calculation from Eq. (^) for a mode with effective index, 
n c ff = (3/k = 1.444. The full and dashed lines show the 
approximations Eqs. (^,^), respectively. 




Fig. 3 

Effective mode-index (solid line, left axis) and effective 
area (dashed line, right axis) of a pcf with d/a = 0.15. 
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Fig. 4 

Half-divergence angle (solid line, left axis) and kw (dashed 

LINE, RIGHT AXIS) OF A PCF WITH d/A = 0.15. 



0.15. The first PCF fabricated by Knight et al. § of 
this kind had a pitch of A = 2.3 fim and was found to 
be single-mode in the range A/A between 1.5 and 6.8. In 
Fig. ^ we show the corresponding half-divergence angle. We 
have also shown the value of the dimcnsionless parameter 
kw (dashed line, right axis); the magnitude justifies the 
application of the approximate result in Eq. (g) to PCFs. 
We note that for non-linear PCFs (!(]] the value of kw will 
approach the regime where deviations from Eq. ([}]) arise. 

In Fig. ^| we show the half-divergence angle for different 
hole sizes where the fiber is endlessly single mode |ll| . For 
small hole sizes d/A we note that in practice the operation 
is limited by a significant confinement loss for long wave- 
lengths where the effective area increases [Q. In Fig. || 
this can be seen as a bending-down of 9 for small A/ A. In 
general the NA increases for increasing hole size and fixed 
pitch and wavelength. By adjusting the pitch A and the 
hole size d this demonstrates a high freedom in designing 
a fiber with a certain NA at a specified wavelength. 

In order to verify our calculations experimentally a PCF 
with d/A ~ 0.53 and A ~ 7.2 /im has been fabricated. In 
Fig. |^ we compare our calculations to a measurement of the 
NA at the wavelength A = 632 nm. As seen the calculation 
agrees well with the measured value. 

IV. Conclusion 

We have studied the numerical aperture (NA) of pho- 
tonic crystal fibers (PCF). The calculations is based on the 
approximate "standard" result 9 ~ tan _1 (A/7rw) which we 
have found to be valid in the regime relevant to PCFs. As 
an example we have applied it to the fiber first fabricated 
by Knight et al. j|] . By studying the effect of different hole 
sizes we have demonstrated that the PCF technology have 
a strong potential for NA-engineering in the single-mode 
regime. 



Fig. 5 

Half-divergence angle of a PCF for different hole sizes. 
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Fig. 6 

Half-divergence angle of a PCF with d/A ~ 0.53 and 
A ~ 7.2 /mi. The solid line is a calculation based on the ideal 

STRUCTURE AND THE DATA POINTS ARE MEASUREMENTS AT 
A = 632 nm, 780 nm, AND 980 nm OF THE FIBER SHOWN IN THE 
LOWER RIGHT INSERT. THE UPPER LEFT INSERT SHOWS A NEAR-FIELD 
IMAGE AT A = 632 nm. 
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